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1. INTRODUCTION 


Recently, the heterogeneous catalytic oxidation of 
arenes namely, cumene and tetralin and the cyclic 

comoound, cyclohexene, have been subjects of numerous investi- 
gations [l-lO]. The mechanisms of these reactions have drawn 
considerable academic and industrial attention. In the past 
five years there has been consideraole advancement in the 
elucidation of the mechanisms of oxidations of cumene and 
cyclohexene. A brief literature survey on the mechanisms of 
oxidation of these compounds is presented in the following 
section: 

1 , 1 Literature Survey ; 

Mukherjec and Uraydon [1] have studied the liquid- 
phase oxidation of tetralin (1,2,3,4-tetrahydronaphthalene) 
using both soluble and insoluble catalysts tuid compared the 
reaction rates. Oxides of nickel, manganese, and copper were 
found to be extremely active while those of aluminium and 
zinc were inactive. Rates of oxidation were measured in the 
temperature range of 45-90°C. The initial product of oxidation 
reaction was found to be tetralin hydroperoxide which decomposed 
further into tetralone and tetralol. A ketone: alcohol ratio 
of 2:1 was observed in the product with most of the catalysts. 

It was further noticed that a critical hydroperoxide to 
catalyst ratio existed below which the reaction did not proceed. 
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rlin.'tic studies have been made with four of the best catalysts, 
nn:i tne rrnction mechanism proposed is as follows: 


Ini tie t ion 

ki 


ROOH’M — 
RO‘M + RH 
•OliM + RH 



RO*M + *0HM 
ROH + R- 
H 2 O + R* 


Pro pu. Ration 

k . 

kp 

HO^ + RH —2 — > 


HO- 

aooE 


+ R* 


T(-rmirntion (bimolecular) 
k 

2 HO* ^ > 2R" CO + H 2 O + O 2 

Termination (monomolecular ) 

k 

RO* 1— > R” CO + 0^ + H.O 

^ k 2 2 

RO^ 0 ^ J^I. QQ 

wi til R C^.H^C(CH^)2 and M -- active site of the 

catalyst surface 

Van iiam et al. [2] have found that the use of silver 
catalyst siiowud a significant lowering of the activation energy 
in thf' oxidation of cumene. They also found a pronounced 
increase in the oxidation rate by using silver instead of 
silvt r~gold alloy as the catalyst. 

Hf.-'uburg et al. [3] studied the heterogeneous liquid- 
phastf oxidation of cyclohexene using manganese dioxide as the 
catalyst. The decomposition of cyclohexenyl hydroperoxide was 
also studied. The dependence of initial rate on the catalyst 
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'■(nount h'li hx‘r<n studied both for oxidation and decomposition 
re‘’.ctionR r.nd def^r-neratc- chain branching mechanism had 
b» *. n rostul’ited by the authors. Hydroperoxide which was 
produced in the propagation step was decomposed on the solid 
generatinf the free radicals necessary for the initiation of 
t;.c chains. Cyclohexonol and cyclohexenone were produced in 
tht initiation and termination steps. Above a critical catalyst 
to hjairocarbon ratio, no oxidation reaction was observed. At 
th-'t critical composition, the rate of hydroperoxide decomposition 
was i;que.l to its rate of formation. The length of the chain 
d.,:cr<:*ased to a low value with increasing catalyst to hydrocarbon 
ra.tios. Th-.- mechanism adapted is as follows: 

Initia.tion 


ROOH + M - 
KH 4 K..H0 
RO* + RH - 



RO* + M..HO* 
■™» rr+M+H^o 
ROH + R* 


Pron'‘.,-:ation 


k. 

R- + Og ^ RO2 

RO* + RH ^ ROOH + R* 

Homogeneous termination 

k 

RO; + RO; R'=0 + ROH + O 2 

RO* 4- RO* ^ R'=0 + 

The oxidation of cumene and the decomposition of cumene 
hydroperoxide on silver, copper and platinum were studied by 


Casomier ct al. [4]. Cumene hydroperoxide dissolved in dodecane 



r^'iniily it conirtosv d by plritinum and appreciably by silica gel 
C r by silv ,r in tin absence of gaseous oxygen, but virtually 
no d^. cor.rosi tier. w:is observed with unsupported silver in an 
0X5’’ g' n 'ttr.osphorr . The oxidation of cumene to cumene hydro- 
peroxide sV,o.v'i.d a completely different picture; this oxidation 
was strongly catalyzed by silver, whereas copper, gold and 
platinum did not catalyze the reaction. For the cumene 
oxidation v.'i t', silver a free radical (pres\mably ROO' ) was 
dt ti ct< d in the reaction mixture by EPR spectroscopy. The 
r. -ction m^'chanism assumed diatomic oxygen to be adsorbed on 
tie silver surfacf . The proposed mechanism is as follows: 
Initiation 

ROOK RO* +‘0H 

' Or! + HH ^ HgO + R* 

HO* + RH ^ HOH + R* 

ri'opa,‘y tion 

R*02 - ROO* 

HOO* + RH ROOH + R* 

Te-rraination 

2 ROO* ^ non radical product. 

V-irm*'; and Graydon [5] studied the heterogeneous catalytic 
oxid'‘tion of cumene in liquid phase. The investigation was 

“T "T y-T- y 

made- at BO^C with molecular oxygen. Cobalt ' oxide, Mn02, 
NiO and CUgO wore used as catalysts. Detailed studies of the 
('flt'Ct of catalyst weight, hydrocarbon concentration and initial 
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I'oxid' o.'jr.o< ntr.-ition on the rate of oxidation were 

TT TTT 

irti’t'. a •. ;;r,-lo,vin>; cobalt ♦ oxido as the catalyst. The 
cxVf.a r. a.tsorrtion rate was found to be 0.4 power with respect 
to th' c'-. t-’.l.Yst weight, and 1.5 power with respect to the 
onr:.cn‘ cor.Ct ntration. The concentrations of the main products 
cf oxid’-.t ion o"' cuacni. , viz., cumene hydroperoxide, dimethyl 
r!;. nyl carrin.ol ''iid acetophenone wore found to depend on the 
w* i.:ht of t-h- C'’tr*.l.vat. Reaction mechanisms and rate expressions 
:',.r tin oxi:\'-tion of cumene were proposed. The mechanisms are 
illustrated below; 


Initi''tion 
RH + M 



R* 


Fropa*'”. tion 
H* -i 0., 
HH 



ROO* 

ROOH + R* 


?' rm 2 nation 

2i00’ + ROO' 


2i00* + ROO* inactive products 

7ht 'I't 0 V(.;‘ mechanism was postulated when there was no 
hyironc roxid.- present in cumene. The effect of the initial 
hy ircmeroxidi, would be felt in the initiation reaction and a 
isoinej bli- raecl.anism proposed wasi 


R00H....M — 
RO* M + RH 


RO* M + *0H M 

ROH + R* 


'OH M + RH 


H2O + R* 
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Neuburg et al. [6,7] conducted a detailed investigation 
on th.e oxidation of cyclohexene and decomposition of cyclo- 
hexene hydroperoxide. The liquid -phase decomposition as also 
the oxidation of cyclohexen,e? hydroperoxide in cyclohexene [6] 
had been studied in the temperature range of 40 - 70*^0 using 
three typos of Mn 02 with different specific surfaces as 
catalysts. The influence of initial hydroperoxide concentration 
catalyst weight to liquid volume ratio and temperature were 
studied in each case. It was postulated tha.t the hydroperoxide 
formed an equilibrium complex with the catalyst which further 
decomposed to yield free radicals. The initial activity of 
oxides decreased during the reaction and the observed decay 
in the rate of hydroperoxide decomposition ha.d been attributed 
mainly to the deactivation of catalyst surface. The decompo- 
sition products wore cyclohcxenol, cyclohexenono and water. 

The number of molos of alcohol formed was higher than the 
number of molos of ketone in each case. In the oxidation 
study, reaction periods upto 60 min. were used to study the 
initial activity of the cata.lyst. The maximiun rate of cyclo- 
hexene oxidation was 0.5 order with respect to the catalyst 
in the region of low catalyst weight to hydrocarbon ratios. 

A much lower order was found for higher ratios. Above a 
critical amount of catalyst, the oxidation was comnletely 
inhibited due to the quantitative decomposition of hydroperoxide. 
The maximum oxidation rate was proved to be first order with 
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on the oxidation of cyclohexene and decomposition of cyclo- 
hexene hydroperoxide. The liquid-phase decomposition as also 
the oxidation of cyclohexen.e? hydroperoxide in cyclohexene [6] 
had been studied in the temperature range of 40 - 70*^0 using 
three types of Mn 02 with different specific surfaces as 
catalysts. The influence of initial hydroperoxide concentration 
catalyst weight to liquid volume ratio and temperature were 
studied in each case. It was postulated that the hydroperoxide 
formed an equilibrium complex with the catalyst which further 
decomposed to yield free radicals. The initial activity of 
oxides decreased during the rea.ction and the observed decay 
in the rate of hydroperoxide decomposition had been attributed 
mainly to the deactivation of catalyst surface. The decompo- 
sition products wore cyclohoxenol, cyclohexonono and water. 

The number of moles of alcohol formed was higher than the 
number of moles of ketone in each case. In the oxidation 
study, reaction periods upto 60 min. were used to study the 
initial activity of the catalyst. The maximum rate of cyclo- 
hexene oxidation was 0.5 order with respect to the catalyst 
in the region of low catalyst weight to hydrocarbon ratio s- 
A much lower order was found for higher ratios. Above a 
critical amount of catalyst, the oxidation was comnletely 
inhibited due to the quantitative decomposition of hydroperoxide. 
The maximum oxidation rate was proved to be first order with 
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rosTDCct to tho concentr-T.tion of cyclohexene. The oxidation 
rate? x^ras of zero order with resrect to the initial hydropero- 
xide concentration in the range of 3 x 10 ~^ to 0 . 5 gmole/liter , 
giving support to the theory that free radicals were produced 
from the slov; decomposition of the catalyst-hydroperoxide 
complex. Oxidation products such as cyclohexenol and cyclo- 
hext-'none innibited the oxidation rate when they were initially 
nresc-nt in concentrations greater than 50 pmole/ml. Apparent 
'ictivation onorgios for the oxidation ranged between 11.4 
and 13.0 Kc-xl/molc. Product distribution studies showed that 
cyclohexenol '.'.nd cycloxenone were produced in nearly equal 
amounts. A degenerating chain branching mechanism was proposed 
in which the catalyst played an important role in the initiation 
and tfrmin'' tion stops. Rate equations derived from this 
mt on", tii cm ox l'iin<.'d quantitatively the experimental observations. 

I’anayo lova ot al. [ 9 ] worked on the role of metallic 
enpnor in thi* oxidation of cumene to cumene hydroperoxide. 

They reported that at moderate temperatures copper protects 
cument! hydronaroxide against decomposition. Copper corrosion 
takes plr’ce upon continuous processing of raw material, parti- 
cul'xrly .-It high temperatures. Interaction between the metal 
extracted into reaction modixim and cemiene hydroperoxide lead 
to intemso decomposition of cumene hydroperoxide. 

Srivastava and Srivastava [ 8 ] investigation the liquid 
phase oxidation of cumene with Cr^O^, Mn 02 and ^'© 2*^3 catalysts 
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in the temperature range of 70 - 90°C.’ Similar to the earlier 
finding of Varma a.nd G-raydon, the present authors observed 
a critica.1- hydroperoxide to catalyst ratio, below which the 
reaction did not proceed. Effects of catalyst weight and 
hydrocarbon concentration on the rate of oxidation were studied 
in detail. The apparent activation energy for the overall 
oxidation was calculated to be 9.8 Kcal/g mole. A reaction 
mechanism had been proposed and the rate expression derived 
was in good agreement with the experimental data. 

Agarwal and Srivastava [10] studied the kinetics and 
catalytic behaviour of NiO supported on alumina and activated 
carbon and also Ni MoO^ supported on alumina in the liquid 
phase oxidation of cumene. The catalysts had been characterized 
by means of DTA, TGA, X-ray, IR and BET studies. Detailed 
studies of the concentrations of main products along with 
off6:!ct of catalyst weight, and hydrocarbon concentration on 
the rate of oxidation were conducted at 80^0 employing these 
supported catalysts. A re'^'ction mechanism was proposed in which 
the catalyst plays an important role in the initiation and 
termination steps. Rate equations derived were in good agree- 
ment with the experimental data. It was observed that use of 
carbon as a support instead of alumina (NiO catalysts) induced 
increased activity, lower activation energies, and shifted the 
limiting oxidation rate to lower catalyst amount. 
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It has been established that, for the heterogeneous 
liquid phase oxidation of cumene, cyclo hexene and tetralin 
catalyzed by transition metal oxides, there is a critical 
weight of catalyst above which the rate of oxidation drops 
ca tastrophicalljj- [1,3, 5,8], which suggested that the rates of 
formation and destruction of hydroperoxide became equal [7]. 

A general mecha.nism explaining the simultaneous formation and 
destruction of hydroperoxide in the oxidation of cyclohexene 
has recently been proposed by leuburg et al. [7], and used 
by Agarwal and Srivastava [10] to explain the kinetics of 
cumene oxidation. However, very few data have been presented 
for these systems in the region where inhibition was effective. 

1.2 Present ¥ork : 

In tho present work MnO^ catalyzed oxidations of cumene, 
cyclohoxenc and tetralin in tho region where inhibition was 
effective, have been studied at the temperatures of 80, 60, 
and 65°C respectively. The results favor the view that these 
catalyzed oxidations involve simultaneous formation and destructioj 
of hydroperoxide. Also the chain branched character of the 
generation of free radicals is confirmed. 



2. EXPBRIMMTAI 


2.1 Materials ; 

Cumene, cyclohexene and tetralin were obtained from 
Herdillia Chemicals, Koch-light (England) and Riedel De Haen 
Ag (Germany), respectively. Cumene and cyclohexene were 
further distilled fractionally. 

Monochlorobenzone of Sarabhai M. Chemicals was used 
as inert solvent for varying the concentration of hydrocarbons 

Hydroporoxides of cumene, cyclohexene and tetralin 
were prepared separately by the thermal oxidation of the 
respective hydrocarbons at 60°G using MnO^ as the catalyst. 

Oxygen and nitrogen were supplied by Indian Oxygen 
Limited, Kanpur. Nitrogen was used as an inert for the decompo 
sition runs. 

MnO^ catalyst was supplied by BDH (India) Limited, 

Bombay. 


2 , 2 Apparatus ; 

The apparatus w.as essentially similar to that used by 
Bolland [ll] for the thermal oxidation of ethyl linoleate. 

The apparatus consisted of a graduated gas-burette. 
The burette was 50 cm long and 1 cm inner diameter with both 
ends open. The bottom was connected to a mercury reservoir 
through a polythene tubing as shown in Figure 1. 
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KEY TO FIGIXRS 1 

B 55/54 Glass Joinl Reactor 
Constant Temperature ¥ater Bath 
Remi Stirrer 
Rubber Tubing 
Three-way Stop Cock 

Connection to Vacuum Rump through Trap 

Mercury Manometer 

Gas-Burette With Mercury 

Polythene Tubing 

Mercury Reservoir 

0:^alic Acid Container 

12V Battery 

Connection to Gas Cylinder 
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Oxygen line and the reactor were connected to the 
gas-burette by polythene tubing through a three-way stop- 
cock. The reactor consisted of a B 55/54 glass joint. The 
male joint was connected to a three-wa3/- stop-cock, which in 
turn was connected to a va.cuum lino and mercury manometer. 

The reactor volume was 400 cc. A vacuum pump was used for 
evacuating tho system. The reactor was placed in a constant 
temperature water bath. ¥ater in tho bath was stirred by a 
^emi-stirrer . The bath was heated by a 2 K¥ heater. An on- 
off relay and a, contact thermometer of range 0-500^0 were 
used to maintain the desired bath temperature. 

Efficient stirring of the reaction mixture was achieved 
by using a magnetic stirrer. The stirrer consisted of a 
glass-covered iron piece. 

2.5 Procedure ; 

All the oxid'^tion runs were performed using MnO^ as 
the catalyst at the corresponding temperatures mentioned 
below for different systems. 

Cumene 80°C 

Cyclohexene 60°G 
Tetralin 65^C 

The experimental procedure consisted of the following 

steps: 
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1. The water-bath was switched on and allowed to reach 
the desired steady-state temperature, which normally took 
half an hour. 

2. The desired catalyst amount was taken in the clean 
and dry reactor. Then the iron stirring piece was dropped in 
the reactor. 

3. 25 ral of the hydrocarbon were added in the reactor. 

4. The reactor was connected to the gas-burette and 
the system was evacuated by means of a vacuum pump. 

5. The system was filled with oxygen. 

6. The reactor was placed in the bath. 

7. After fifteen minutes the ma,gnetic stirrer was 
switched on and change in the burette reading was recorded. 

8. After the run the products wore allowed to cool 
down to room temperature and then centrifuged for fifteen 
minutes at 800 RPM in order to separate the catalyst from the 
products. 

Hydroneroxido Decomnosi tion Runs : 

For the decomposition runs, the hydroperoxide solution 
of the respective hydrocarbon of known concentration was 
added, together with the desired amount of MnO^ catalyst to 
the reactor, under nitrogen. All the decomposition runs were 
performed for thirty minutes. 

The samples were analyzed as described below: 



2 . 4 Aaal.ytical 


The main products of decomposition were the respective 
hydroperoxide, alcohol and ketone. In the case of cumene, cumene 
hydroperoxide, cumyl alcohol and acetophenone were the products. 
Similarly for the decomposition of cyclohexene, cyclohexene 
hydroperoxide, cyclohexenol aixa cyclohexenone were the products, 
while tetralin hydroperoxide, tetralol and tetralone were the 
products in the decomposition of tetralin. The following methods 
were used for analyzing the various products of the decomposition 
runs. 

hydroperoxide : 

The hydroperoxide content in each ca,se was estimated by 
iodometric method [12], 40 ml of isopropyl alcohol and 5 ml 

sample were mixed and refluxed. The mixture was cooled to room 
temperature and 10 ml isopropyl a,lcohol saturated with sodium 
iodide were added aid refluxed. The liberated iodide was 
titrated against sodium thiosiibphat e solution. 

Alcohol; 

Cmnyl alcohol aid peroxide were converted to phenol by 
treating the sample with sulfuric acid and hydrogen peroxide[l3] , 
Bromide-br ornate method was used to estimate phenol [17] . 5 ml 
of the sample was treated with 50 ml of O.IB of potassiitm 
bromat e-bromide solution and 5 ml of HOdi*. The mixture was 
shaken and allowed to staid for 15 min. Two grams of potassium 
iodide were added. The liberated iodine was titrated with 0.1¥ 
sodium thiosulphate using 2 ml of starch indicator. 

Ketone ; 

The ketone concentration in the product was detenained 
using Berkin-Elmer- 137 IB. spectrophotometer by studA^ing the 
absorption at 1700 cm"”^. 



3. RESULTS AMD PI SOUS SI ON 


The symbols RH, ROOH and ROH are used in general for 
hydrocarbons, hydroperoxides and alcohols, respectively. All 
the relevant experimental data for the oxidation and decompo- 
sition runs are presented in Appendices A, B, and G respectively. 

In order to study the effect of hydrocarbon concentration 
in the region where in|iibition was effective, it was necessary 
to study the oxidation rate as a function of the catalyst 
ratio for all the three systems. The catalyst ratio is defined 
as the weight of the catalyst in grams divided by the volume 
of the hydrocarbon in ml. 

3 . 1 Initiation s 

The reaction did not commence when the cumene oxidation 
was performed at 80°C in absence of the catalyst or any free r 
radical initiator. When some hydroperoxide was added to the 
system the oxygen consumption started immediately. Therefore, 
to initiate the reaction both catalyst and hydroperoxide were 
necessary in all cases. Hydroperoxide concentration to catalyst 
ratio of approximately 2.5 x 10 ^ graole/lit g was used in all 
the three cases. 

3 . 2 Effect of Catalyst Weight : 

Oxidation of Cumene : 

The amount of oxygen absorbed versus time of oxidation 
for cumene at 80°C is shown in Figure 2. Data were obtained 
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on the rate of oxygen a,bsorption for various catalyst weight 
to hydrocarbon volume ratios ranging from 0.02 to 0.2 gram of 
catalyst per ml. of hydrocarbon. The initial rate obtained 
for 0.02 g/ml was 1.7 x 10~ g.mole/sec. Srivastava and 

—7 

Srivastava [8] observed an oxygen absorption rate of 1.3 x 10 ' 
g mole/sec under the same conditions. 

The rate of oxygen absorption versus catalyst weight 
is shown in Figure 3. It was observed that the rate of oxygen 
absorption was maximum at 0.14 g/ml. A further increase in 
the catalyst weight resulted in the decrease of rate of oxygen 
absorbed. This behaviogur was observed by Srivastava and 
Srivastava [8] as well as by Varma and G-raydon [5]. The value 
above which the rate declines is referred to as the critical 
catalyst amount. Srivastava and Srivastava [8] established 
that the critical catalyst amount was unaffected by variation 
in tempera. ture over the temperature range 70 - 80°0. 

Oxido-tion of Cyclohexene : 

The data recorded in Table B.l of Appendix B represent 
the consumption of oxygen with time for various catalyst ratios 
with Mn02 as a cntalyst at 60°C. A detailed study was performed 
to obtain the rate of oxygen absorbed for various catalyst 
ratios ranging from 0.08 to 0.4 g/ml. Figure 4 shows the 
rate against catalyst ratio. As in the case of cumene the rate 
of oxidetion dropped rapidly after 0,2 g/ml. Similar 
observation was made by Beuburg et al. [3]. They observed a 
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decline in r-ite at 0.16 g/ml under the same conditions. In 
the present study the maximum rate obserred was 3 x lO"'^ g.mole/ 

— *7 

sec. while Neuburg et al. [3] found 7 x 10~ g. mole/sec. 
Oxidation of Tetralin ; 

The oxidation of tetralin was carried out using the 
same Mn02 ca,talyst at 65°C. The amount of oxygen absorbed 
versus time is shown in Table G-1 (Appendix O). Oxygen absor- 
ption was studied at catalyst ratios of 0.2, 0.4, 0.8, 0.9 and 
1.0 g/ml. The rate data were plotted against catalyst ratios 
as shown in Figure 5. As seen the behaviour is similar to 
those of oxidation of cumene and cyclohexene. Mukherjee and 
Graydon [l] also observed the same behaviour while investigating 
this reaction. 

Thus, in all the three cases of oxidation of cumene, 
cyclohfixeno, and tetralin there is a catalyst ratio above 
which oxidation rates decline. It has already been established 
that before the break -point the rates of oxidation in the 
presence of transition metal oxides proceed with anproximately 
0.5 and 1 power with respect to catalyst ratio and hydrocarbon 
concentration, respectively. [2,5,8,10] . 

5 • 3 Effect of Hydrocarbon Cone entra.t ion in the Inhibiting 
Region : 

The influence of cumene, cyclohexene and tetralin 
concentrations on the respective rate of oxidation was studied 
with the mixture of the corresponding hydrocarbon with 
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monoclilorobenzene at catalyst ratios chosen from the inhihitir, 
region. 

Cumene ; 

Runs were conducted for three different catalyst ratios 
0.16, 0.18 and 0.2 g/ml (Figure 3). Figure 6 represents a 
typical plot of oxygen absorbed versus time at various hydro- 
carbon concentrations for the catalyst ratio 0.16 g/ml. The 
data for 0.18 and 0.2 g/ml. are presented in Table A-2 (Append 
A) . The rates were calculated and plotted against the cmaene 
concentration, on a log-log graph as shown in Figure 7. It 
is interesting to note that the orders obtained were 2.3, 2.1, 
and 2.2 with respect to hydrocarbon concentration for the 
catalyst ratios 0.16, 0.18 and 0.2 g/ml, respectively. It 
shows that the rate dependence with respect to the cumene 
concentration in the inhibiting region is approximately 2. 

Cyclohexeno ; 

As in the case of cumene the cyclohexene concentration 
was varied and runs were conducted for four different concen- 
trations. The catalyst ratio chosen was 0.22 g/ml (Figure 4). 
Table B.2 (Appendix B) represents the data of oxygen a.bsorbed 
versus time. The rate of oxidation as a function of cyclo- 
hexene concentration is shown in Figure 8 on a log-log plot. 
The rate dependence with respect to cyclohexene concentration 
was found to be 1.7. 




rves for catalyst ratio 0.15 g Mn02 pe 
tction of cumene concentrations qmoU 
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Trtrnlin : 

In caso of totralin a ca.talyst ratio of 0.9 g/ml was 
employed. This was again chosen from the inhibiting region 
(Figure 5). The oxygen absorption data, are shown in Table C.2 
(Apnendix C). In Figure 9 rates of oxygen absorption are 
plotted against hydrocarbon concentrations. on a log-log plot. 

The order with respect to hydrocarbon concentration was found 
to be 2.1. 

3 . 4 Decomposition of Hydroperoxide ; 

The decomposition reactions of the hydroperoxides of 
cumene, cyclohoxene and tetralin were studied for various 
catalyst ratios at the temperatures of 80, 60 and 65*^0, 
respectively. The decomposition products observed were alcohol, 
kotonc and water. The data are presented in Tables A. 3, B.3 
and C.3. Figures 10-12 represent the percentage 
of hydroperoxide decompos;fe‘d ■ against various catalyst ratios 
using a single initial hydroperoxide concentration in each case. 
It is interesting to observe that the percentage decompostion 
increases with the increasing catalyst ratio upto the break-in- 
point and then shows a decrease on further increase in the 
catalyst ratio. 

3.5 Reaction Mechanism ; 

The experimental observations made during the study 
of oxidation of cumene, cyclohexene and tetralin may be 








nitial hydroperoxide =44 jU moie/T’t 
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summf'.rized in general as follows: 

(a) For lovjer catalyst ratios the rate of oxidation 
increa,ses with increasing catalyst ratio. For 
higher catalyst ratios, there is a decline in 
the rate of oxidation which is referred to as the 
inhibiting region, 

(h) For the catalyst ratios in the inhibiting region, 
the rate of oxidation is a.pproximately second 
order with respect to the hydroca.rbon concentration. 

(c ) The decomposition of hydroperoxide in all the 
three cases increases with increasing catalyst 
ratio, passes through a maximum and then declines 
suddenly. 

Based on the observations made in the present study and 
from thu esto/Dlished facts that, before the limiting region, 
the rate dependence is 0.5 order with respect to catalyst 
weight and first order with respect to hydrocarbon concentration, 
a general mechanism as proposed by Neuburg et al. [7] is 
described below: 

Initiation 

K. 


ROOH + M — 



(ROOH M) 

(1) 


(R00H...M) ^ — > 

R0‘ + M. . .OH" 

(2) 

k. 

ROOH + M. . .0H~ 

k. 

— ^ RO * + H^O + M 

(3) 

RO* + RH ^ ^ 

ROH + R' 

(4) 
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Propagation : 
El* + O 2 
RO* + RH 
Termination; 


RO* + RO* 
RQ* + M 


RO^ 


ROOH + R* 


(5) 

( 6 ) 


k. 


K 


inactive products ^2 


8 


(RO. 


.M) 


k, 


XV. Q 

ROg + (RO 2 ...M) — inactive products 


( 8 ) 

(9) 


3 • 6 Derivation of Kinetic Expression : 

At steady-state, the rate of initiation of free radicals 
is equal to the rate of consiunption of free radicals. 

Therefore; Rate of initiation of free radicals 

R. = 2 k.^ [RO*]^ + kg[R0*][R0*...M] 

Reaction (8) is under equilibrium. 

Therefore I kg = [ggPj ’ 

[ROO M] = kg [RO*] [M] 

R^ = 2 k.^ [RO^]' + kgkg [R0*]2 [M] 


( 10 ) 


( 11 ) 


= [RO*]^ (2 k.^ + kg kg [M]) 


R. 


^ 2k.^ + kgkglM]^ 


1 


( 12 ) 


Rate of oxygen consumption is given by 

- = kg [R0-] [RH] + \ - k^LROj]' 


(13) 


R. 


2 


- ^6^2 k.^ + kgkgtMl^ 

k„ 


+ R.(l 


1 - 2'k.^+kQkg[RH] 


■) 


(14) 
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Ass-oming formation of oxygen through termination is negligible 


d[02] 

dt 


= k 


R. 


6 ^ 2 k^+kQ]£g[M] 


(15) 


If the termination of peroxy radicals on the catalyst 
surface is negligible compared to biradical termination 


^[Op] 

- = ^6 


R. 


2 k, 


7 


1 


[RH] 


( 16 ) 


It may be assumed that R^ is proportional to the amount 
of catalyst. This rate expression (Eq. 16) clearly accounts 
for the observations of 0.5 and 1 order with respect to the 
catalyst ratio and the hydrocarbon concentration, respectively. 

When the limiting rate condition is reached, the rate 
of formation of hydroperoxide by reaction (6) becomes equal to 
the rate of decomposition on the catalyst surface, R^. Hence 
the net rate of formation of hydroperoxide is 


let 


d[R00H] 

dt 

R. = a R 


~ ^6 ^2 k^+kQkg[M']^^ - - 0 


1 - -^d 

where, a is the fraction of hydroperoxide decomposed. 

From equations (17) and (18) , solving for R^ give 


(17) 

(18) 


k" a" [RH]^ 


\ " 2k^ + k^kjl] 


(19) 


^8 9>- 


By substituting the value of R^^ in Equation (14 ) , the 
limiting rate of oxygen consumption becomes 
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d[02] cc k| [RH]^ 

= ____^] [i+a(i - 2k^+kQkg[M]^ ^ 

( 20 ) 


For small value of a the limiting rate of oxidation is given "by 


d Op] 

_ ( 2_n 

^ dt ^co 


a kg [RH]^ 

2k^+kgkg[M] 


( 21 ) 


Therefore, under limiting condition the oxidation rate 
must show a second order dependence on the hydrocarbon concen- 
tration. The experimental observations for all the three 
systems cumene, cyclohexene and tetralin give a rate dependence 
of approximately 2. The observations are, therefore, in good 
conformity with the rate expression derived. The similarity 
in the kinetic behaviour of the three systems may be dixe to 
formation of resonance stabilized radical species. Therefore, 
it is envisaged that a single mechanism satisfactorily explains 
the kinetic behaviour observed for all these systems. 



4. STATISTICAL ANALYSIS OF THE DATA 


In order to further substantiate the model a statistical 
analysis was done. The statistical treatment involved esti- 
mation of parameters for the three systems under study, 
establishing confidence intervals [14] for the parameters and 
then a residual analysis [15]. 

4.1 Estimation of Parameters and Their Confidence Intervals : 

The parameter estimation was essentially based on the 
minimization of the sum of squares between the expected and 
the experimental responses with a suitable convergence criterion. 
Rosenbrock [16] optimization program was used for the minimi- 
zation. IBM 1800 and IBM 7044 computers were used. The method 
is of direct-search type. An initial estimate was used to 
start with. The program was tested for convergence using 
Banana function, 

E = 25 (1-x^)^ 

For a fast convergence and simplification the non-linear 
model was reduced to a three parameter one as follows; 

^[0^] b^ [RH]2 
“ dt = b2+b^[M] 

where, b^ = a kg 

bg = 2 

bj = kg kg 




37 


The problem was reduced to estimating the three 

parameters b^, b^ and b^ in each of the three systems. 

The parameters thus obtained and their 95 per cent 
confidence intervals for the three systems are tabulated below. 
The program listings and the calculation of confidence intervals 
are included in Appendix Di.. 


System 

Parameter Values 

95 per cent confidence 
intervals 


Cumene 

bi 

= 

2.2588 

X 

10^ 

1.9598xlO,^<: 

p^c. 2.5578 

X 

10^ 


h 


1.3744 

X 

10^° 

1.1924xl0^i 

1.5564 

X 

lo^o 



= 

1.2048 

X 

10^ 

-1.0600x10^?^ 

p^< 1.0600 

X 

IqIO 

Cyclohexene 

h 

= 

2.1638 

X 

10^ 

1.7758xlO^:fc 

2.5518 

X 

10^ 


h 

= 

1.3744 

X 

10^0 

1.1274x10^£ 

p2< 2.6214 

X 

lolo 


’’3 


1.0500 

X 

10^ 

-1.1280xl0^:§ 

P^v- 1.1200 

X 

IqIO 

Tetralin 

h 

= 

1.9584 

X 

10^ 

-2. 1042x10^ <- 

p^< 2.5058 

X 

10^ 


h 


3. 6280 

X 

10^ 

3.0140xl0^e 

p2<v 4.2419 

X 

10^ 


h 


1.8140 

X 

10^ 

I,l326xl0^4i 

p^< 2.4954 

X 

10^ 


An examination of the confidence intervals indicates 
that the parameter values are all positive in the confidence 
region. This is one of the importa.nt criteria for model 
selection. 

4 . 2 Residual Analysis : 

The residual analysis is a simple but powerful tool 
to test the adc?quacy of a model. A residual is the difference 
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"between an o"bserved value and the corresponding predicted 
value, "Jsing the estimated parameters, rates were predicted 
(Appendix D3) . The residuals were then obtained by subtracting 
the predicted rates from the observed ones. Figures 13, 14 
and 15 show the residual plots a.gainst the independent variable 
viz. hydrocarbon concentration for the systems cumene, cyclo- 
hexene and tetrrolin, respectively. The plots show a random 
behaviour of the residuals which further substantiates the 
model. 



Residual x 10 



Fig. 13 'Plot of the r esrluals ■ ugain^-t the hydrocat bo^ 

concentrotion for the system cumene. 




















5. CONCLUSIOIS 


The kinetics of oxidation of the three systems - 
cumene, cyclohexene end tetralin showed a similar behaviour. 

In all the three cases it was established that as the catalyst 
to hydrocarbon ratio wa.s increased, the rate of oxygen 
absorption increased and after a certain ratio the rate 
declined suddenly. While an order of one with respect to 
hydrocarbon concentration was established earlier in the 
region before the break-in-point, a second order rate dependence 
on hydrocarbon concentration was obtained in the inhibiting 
region. The second order dependence in the inhibiting region 
supports the degenerate chain-branching mechanism. 

The observed data were given a statistical treatment 
to test the adequacy of the model. The parameters were 
estimated using Rosenbrock's optimization method. The confi- 
dence intervals were also estimated. The parameters were all 
positive in the confidence region. Finally, a residual analysis 
showed a random nature with respect to the independent variable, 
hydrocarbon concentration. The positive nature of the parameters 
and the random behaviour of the residuals substantiate the 
model derived on the basis of degenerate chain-branching 
mechanism. Thus, such model, in general, accounts for the i 

reactions where hydroperoxide forms and propagates the chain 
reaction. 
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APPENDIX 


OOTIENE RUNS 



q. D 


TABLE A-1*. OXYGEK ABSORPTION VERSUS TIME AT DIFFERENT MnO 

CATALYST WEIGHTS AT 80°C 


Run No . 


Catalyst 
ratio , 
g/ml 


Time , 
min. 


Oxygen consumption 
(g mole/ml)xlO^ 


Initial rate, 
g mole/sec 

X 10^ 


C-1 0.02 0 0 


5 

0.602 

10 

1.160 

15 

1.722 

20 

2.250 

25 

2.940 

30 

3.364 

35 

3.785 

40 

4.080 

0 

0.000 

5 

1.080 

10 

2.065 

15 

2.940 

20 

3.505 

25 

3.880 

30 

4.020 

35 

4.235 

40 

4.315 


1.934 


3. 600 





Table A-1 (contd) 


1 


Run No . Catalyst Time , 

ratio, min. 

g/ml 


Oxygen consumption, Initial rate 

(g nole/ml)x 10^ xio'^ 

s ©c 


C -3 


C -4 


C -5 


0.10 0 0.000 

5 1.340 

10 2.412 

15 3.205 

20 3.719 

25 4.100 

30 4.290 

35 4.480 

0.14 0 0.000 

5 1.763 

-10 3.05 5 

15 3.840 

20 4.310 

25 4.520 

30 4.660 

0.160 0 0.000 

5 1.384 

10 2.550 

16.5 3.449 

20 3.719 

25 4.025 

.30 4.220 

35 4.340 

40 4.440 


4.170 


5.170 


4.800 
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Table A-1 (contd) 


Run No. Catalyst 
ratio , 
g/ml 


Time , 
rain. 


Oxygen consumption, Initial rate, 
(g mole/ml) xlO^ g mole ^ ^q7 


C-6 0.18 0 

5 

10 

15 

20 

25 

C-7 0.20 0 

5 

10 

15 

20 

25 


0.000 4.580 

1.276 
2.389 
3.110 
3.460 
3.745 

0.000 2.866 

0.814 

1.850 

2.429 

2,879 

3.140 
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TABLE A-2; OXYO EK ABSORPTION VERSUS TUO FOR DIFFERENT 
CUMENE CONCENTRATIONS (LIMITING ZONE) 


Run Catalyst 
No. ratio, 
g/ml 


Cumene 

concen- 

tration, 

g mole/lit 


Time, 

min. 


Oxygon absor- Initialrate, 
P'tion, g mole ^ ^q7 

g mol e ^ ^q 4 sec 
ml 


C-8 0.16 


G-9 


C-10 


7.82 0 0.000 

5 1,384 

10 2.550 

16.5 3.449 

20 3 . 719 

25 4.025 

30 4.220 

35 4.340 

40 4 . 410 

6.26 0 O.COO 

5 1.078 

10.5 2.015 

17 2,905 

20 3.175 

25 3-519 

30 3.685 

35 3.920 

4.69 0 0.000 

5 0.414 

10 0.909 

15 1.528 


4.800 


4.130 


1.731 
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Table A-2 (contd) 


Run 

i\: O . 


Catalyst 

ratio, 

g/ml 


Cumene con- 
centration, 

g mole/lit 


Time, 

min. 

Oxygen 
ption, 
g mole , 

ml 

20 

1.928 

25 

2.514 

30 

2.930 

35 

3.392 

40 

3.760 

0 

0.000 

5 

0.299 

10 

0.644 

15 

1.018 

20 

1.390 

25 

1.738 

30 

2.049 

35 

2.375 

40 

2.680 

0 

0.000 

5 

0.197 

10 

0.380 

15 

0.583 

20 

0.814 

25 

1.031 

30 

1.261 

35 

1.471 

40 

1.710 


Initial 
rate , 

g mole^ ^j^Q? 


sec 


C-11 


3.99 


C-12 


0.16 


3.13 


1.142 


0.667 



Table A-2 (contd) 


Run Catalyst Cument-^ con- 
No. ratio, centration, 

g/ml g molo/lit 


0-13 0.18 7.82 


G-14 


6.26 


C-15 


4.69 


Time, 

min. 

Oxygen absor- 
ption, 

« mole.. t« 4 
ml ^ 

Initial 

rate, 

« mole^ 

sec 

0 

0.000 

4.380 

5 

1.276 


10 

2.389 


15 

3.110 


20 

3.460 


25 

3.745 


0 

0.000 

3.430 

5 

1.004 


10 

1.940 


15 

2. 640 


20 

3.230 


25 

3.550 


30 

3.810 


35 

3,990 


40 

4.120 


0 

0.000 

2.166 

5 

0.634 


10 

1.229 


15 

1.754 


20 

2.209 


25 

2.691 


30 

3.135< ^ • 

>' t 1, h' 

35 


iBRAftY 


Att. N®. A 
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Table A-2 (contd) 


Run 

lo. 

Oatalyst 
ratio , 
g/ml 

Cumene con- 
centration, 
g mole/lit 

Time, 

min. 

Oxygen ab- 
sorption, 
g mole .q 4 
ml 

Initial 
rate , 

S-mole 10? 

sec 

0-16 

0.18 

3. 99 

0 

0.000 

1.220 




5 

0.345 





10 

0.732 





15 

1.173 





20 

1.532 





25 

1.864 





30 

2.154 





35 

2.470 





40 

2.749 


0-17 


3.13 

0 

0.000 

0.667 




5 

0.207 





10 

0.394 





15 

0.593 





20 

0.828 





25 

1.008 





30 

1.256 





35 

1.470 





40 

1.711 


0-18 

0.20 

6.26 

0 

0.000 

3.365 




5 

0.966 





10 

2.070 





15 

2.925 



20 


3.570 
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Table A-2 (contd) 


Run Catalyst Cumene con- m- Oxygen ab- Initial 

No. ratio, centration, ime, sorption, rate, 

g/ml g mole/lit min. g m ole g mol e^^j^Q? 

ml' sec " 


C-19 0.20 


C-20 


4.69 


3.99 


25 

3.990 

30 

4.265 

35 

4.400 

40 

4.500 

0 

0.000 

5 

0.566 

10 

1.160 

20 

2.235 

25 

2.745 

30 

3.130 

35 

3.465 

40 

3.625 

0 

0.000 

5 

0.469 

10 

0.897 

15 

1.241 

20 

1.531 

25 

1.884 

30 

2.090 

35 

2.499 

40 

2.840 


1. 900 


1.500 
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Table A-2 (contd) 


Run 

No. 

Catalyst 
ratio , 
g/ml 

Cumene con- 
centration, 
g mole/lit 

Time, 

min. 

Oxygen ab- 
sorption, 
g mole p,4 

ml 

Initial 

rate, 

g moley_j,Q7 
sec 

0-21 


3.13 

0 

0.000 

0.716 




5 

0.176 





10 

0.401 





15 

0.704 





20 

0.870 





25 

1,090 





30 

1.318 





35 

1.657 





40 

1.815 




TABLE A-3: CUMENE HYDROPEROXI JE DECOMPOSITION AT 80°C EOR 

dieperent Catalyst (MnO„) ratios 
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CyCIOHEXENE RUNS 
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TABLE B-1: OXYGEN CONSUMPTION VERSUS TIME FOR DIFFERENT MnO 
CATALYST WEIGHTS AT 60°C 


Run No. 

Catalyst 
ratio , 
g/ml 

Time , 
min. 

Oxygen absor- 
ption, 
g_mple , ^4 

Initial rate 

sec 

0-30 

0,08 

0 

0.000 

1.335 



5 

0.103 




10 

0.937 




15 

1.756 




20 

2.270 


0-31 

0.12 

0 

0.000 

2,000 



5 

0.571 




8 

0.780 




10 

0.908 




15 

1.274 . 


0-32 

0.18 

0 

0.000 

2.500 



5 

0.586 




10 

0.835 




15 

1.180 




20 

1.370 


1 

O 

0.20 

0 

0.000 

3.000 



5 

0.880 




10 

1.610 




15 

1.990 




20 

2.230 
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Table B-1 (contd) 


Run Fo. Catalyst Time, Oxygen absor- Initial rate, 

ratio, min. ption, g mole ,^^,7 

g/ml g mole „ sec 

ml ^ 


I 

o 

0.32 

0 

0.000 



5 

0.464 



10 

1.215 



15 

1.760 



20 

2.125 

C-35 

0.40 

0 

0.000 



5 

0.461 



10 

1.025 



15 

1.464 



20 

1.859 


2.A00 


le500 
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TABLE B-2- OXYGEN CONSUMPTION VERSUS TIME FOR DIFFERENT GYCLO- 
IIBXBNE CONOENTRATIONS (LIMITING ZONE) AT 60°C 


Run Ho. Catalyst Cyclohexene Time Oxygen Initial 

ratio, concentration, consumption rate 

g/ml g mole /lit g mole .q 4 £_mpl£xlo'^ 

ml sec 


C-36 0.22 


C-37 5.92 


0-38 


7-880 


5.920 


4.930 


0 

0.000 

5 

1.377 

10 

2.254 

15 

2.750 

20 

3.149 

25 

3.380 

30 

3.540 

0 

0.000 

5 

0.615 

10 

1.464 

15 

2.250 

20 

2.855 

0 

0.000 

5 

0.659 

10 

1.518 

15 

1.918 

20 

, 2.430 

25 

2.955 

30 

3.395 


4.90 


2.50 


2.08 
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Table B-2 (contd) 


Run No. Catalyst 
ra-tio , 
g/ml 


Cyclohexene 
cone entnation, 
g mole/lit 


Time, Oxygen con- Initial 
min. sumption, rate 

ml sec 


C-39 


3.945 


0 

0.000 

5 

0.439 

10 

0.922 

15 

1.428 

20 

1.889 

25 

2.320 


30 


2.755 



TABLE B-3; CHYCL OHEUiBE HYDROPERO^aDE DEG0MP03ITI0E AT 60°0 FOR 
DIPEEREET CATALYST (MnO^ ) RATIOS 
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TABLE C-1: OXYGEN CONSUMPTION VERSUS TIME EOR DIPFERENT MnO, 
CATALYST WEIGHTS AT 65^0 


Run 

No. 

Catalyst 
ratio , 
g/ml 

Time , 
min. 

Oxygen consumption, 
g mole „ ^ ^4 
ml ^ 

Initial rate 
S^le ^ 

sec 

C-47 

0.2 

0 

0.000 

3.6 



5 

1.081 




10 

1.750 




15 

2. 300 




20 

2.800 


00 

1 

o 

0.4 

0 

0.000 

6.6 



5 

2.065 




10 

3.500 




15 

4.500 




20 

5.100 

1 

C-49 

0.8 

0 

0.000 

10.1 



5 

1.514 




10 

2. 600 




15 

3. 300 




20 

4. 000 


0 

1 

o 

0.9 

0 

0.000 

9. 6 i 



5 

2.880 

i 



10 

4.110 




15 

5.190 




20 

5.980 


0-51 

1.0 

0 

0.000 

0.600 i 



5 

0.180 




10 

0.450 

! 



16 

1.081 




20 

1.587 



65 


Table C -2 (contd) 

Run Catalyst Tetralin con- Time, Oxygen con- Initial Rai 

Ro. ratio, centration, min. sumption, g mole ^ 

g/ml. g mole/lit g mole, , 

ml 

15 2.129 

20 2.705 

25 3.210 

30 3.675 

35 4.110 

40 4.470 

0-55 3.69 0 0.000 3.00 

5 0.829 

10 1.298 

15 1.730 ; 

20 2.125 

25 2 . f .'25 

30 
35 


40 


2.885 

3.245 

3.605 


O--?- TjiTfiAlUt HYBROEiitOXIDi: DBOOMPOSITI OH AT 65°C BOR DIIBBRJSHT 
CATALYST (MnO,) RATIOS 
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Detected in infrared spectra 
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ST1TISTIC....L AN.;iYSIS 
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I l=I I+l 

PKd I) = (X( 1 )*RH(II)**2) / {2.=frX(2)+X(3)*CAT( I ) ) 

D02CJ=1 ,KN 
CUNTINUE 
FUNXCN=O.C 
n030 1=1, N 

FUNXCf]=FUNXOK+{10.**7*DX(I )-DXl{ I ) )**2 
30 CONTIMUE 
RETURN 
END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 




evaluation of maximln/mininum of a function 8Y rosenbrock method 

V« COMPONENTS OF ORTHOGONAL VECTORS 
SUBPROGRAM FUNXCN DEFIijES FUNCTION F(X,Z) AS F(0,Z) 

X* UNKNOWN VARIABLES 
Z* KNOWN PARAMETERS 
N* NUMBER OF UNKNOWN VARIABLES, X 
NN* NUMBER OF KNOWN PARAMETERS, 2 

NS* NUMBER OF STAGES AFTER WHICH CALCULATION SEQUENCE IS TERMINATED 
MAXMIN = ZERO FOR MAXIMIZATION AND ONE FOR MINIMIZATION 
LIMIT = ZERO FOR FINITE LIMIT AND ONE FOR NO FINITE LIMIT 
ALPHA* A POSITIVE MLLTl PLYING FACTOR FOR INCREASING A SUCCESSFUL STEPROSOOIAO 
BETA* A POSITIVE MULTIPLYING FACTOR FOR DECREASING UNSUCCESSFUL STEPSROSOO l 50 
E* INITIAL STEP SIZE AT THE BEGINNING OF EACH STAGE R0SO0160 

♦ *#=j!*#**#<c^**:;{c*:ic**#X(*:*#*******^**** ROSOOiTO 


ROSOOOlO 
ROS0002U 
ROSOOC30 
ROS00040 
ROS00050 
R0S00060 
ROS00070 
ROS00080 
ROS00090 
ROSOOlOO 
ROSOO 1 10 
ROSOOiZO 
ROSOO 130 


U FL* VALUE OF FINITE LIMIT 

1 F0RMAT(5I5,4F10.6) 

READ 2, (X (I ) ,1=1 ,N) ,( (Vd ,J) , J=i ,N) , 1 = 1 , N) 

2 FORMAT{8F10.61 

PRINT 3, N, NS, E, ALPHA, BETA 


R0S00180 

R0S00190 

R0S00200 

RDS00210 

ROS00220 


3 FORMAT( 1 X,5CH*******************=***************************#***/ 1XR0S00230 

1, lH*',48X,lH*/lXt31H* NUMBER OF VARIABLES IN) *, 16, 12X, IH*/1X , 3R0S00240 
21H* NUMBER OF STAGES (NS) *, IS , 12X , 1H*/1X, 31H* INITIAL TRR0S00250 
3IAL STEP (E) *,F1 .>.6,8X,1H*/IX,12H* ALPHA *, F10.6, 3X, 6HBETA *R0S00260 

4, F 10. 6, 8X, 1H*/1 X,1H*,48X,1H*/1 X,50H’***************************’****ROS00270 



5 jjc ) 

R0S00280 


I FtMAXMIN.EQ.O) GO fO 11 

R0S00290 


I FILIMIT.EQ.O) GO TO 13 

R0S00300 


PRINT 10 

R0S00310 


GO TO 17 

ROS00320 

4 

CALL MINMAX(N,NS,MAXMIN, LI MIT, ALPHA, BETA, E, FL, X, V, Z ) 

R0S00330 


GO TO 8 

R0S00340 

5 

READ 6, (Zd) ,1=1 ,NN) 

R0S00350 

6 

FORMAT! 8F10.6) 

ROS00360 


PRINT 7, (Zd ),I=1,NN) 

ROS00370 

7 

FORMAT! //7X,3HZ ♦,8F15.8) 

RDS00380 


GO TO 4 

ROS00390 

8 

PRINT 9 

ROS00400 
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9 FfiRMAK //47X ,35HCnMf'LETEO REGUIRED NUMBER OF STAGES//) 

STOP 

in FORMAT(1HO,33HMINIM1ZATION WITH NO FINITE LIMIT) 

11 IF (LIMI T.EQ.O) GO TO 15 

PRINT 12 
GO TO 17 

12 FORMAT! 1H0,33HMAXIMIZATI0N WITH NO FINITE LIMIT) 

13 PRINT 14, FL 
GO TO 17 

14 FnRMAT{lHC,35HMIMMIZATI0N WITH FINITE LIMIT FL =,Fi0.6) 

15 PRINT 16, FL 

16 FORMAT! IH0,35HMAXIMIZATI0N WITH FINITE LIMIT FL =,F1D.6) 

17 IF!NN-GT.O) GO TO 5 
GO TO 4 

END 

C 

SUBROUTINE MINMAX ! N , NS , MAXMI N, LI MIT, ALPHA, BET A, E, FL ,X, V, Z ) 
DIMENSION A!4G,40) ,B (40,41) ,C (40,40 ) ,P( 40,40 ),V! 40, 40) 
DIMENSION D(40) ,0(40) ,R (40 ) ,X (40 ) , Y( 40 ) , Z( 40 ) 

S = E 
L = 1 

NSTAGE=-1 

NTRIAL=-1 

NSUXES=0 

NFELUR=0 

no 1 T=X,N 
D! I)=0.0 
D( I) =X{I ) 

R( 1) =0.0 
Y( I) =0.0 
DO 1 J=ltN 
A( I, J)=0. 0 
B(I,J)=0.0 
C{ I, J)=0.0 

1 P(I,J)=0.0 

2 F=FUNXnN(0,Z,N) 

IF{NTRIAL.LT.O) GO 10 9 

IF(MAXMIN.EQ.O) GO TO 15 

I F{LIMIT.EQ. 0) GO TO 18 

3 IF(F.LE.FN) GO TO 13 

4 IF(NSTAGe.LT.O) GO TO 5 

NFELUR=NFELLR+1 
IF(NSUXES.GT.O) GO TO 21 

5 E=-BETA*E 

6 MTRIAL=NTRIAL+1 
Y(L)=E 

DO 8 1=1, N 
0 ( 1 ) = 0.0 
DO 7 J=1,N 

7 0( 1) =0(1 ) + V(I,J)«Y( J) 


R0S0O4IO 

ROS0042D 

RUS004 30 

ROS00440 

ROS0O45O 

ROS00460 

ROS00470 

ROS0O48O 

ROS00490 

ROS00500 

ROGOOSIO 

R OS 00 5 20 

R0S00530 

ROSOOSAU 

R0SO05S0 

ROS00560 

R0S0057D 

ROS0D580 

ROS0059Q 

R0S00600 

R0S006iU 

R0SC0620 

ROS00630 

RDS00640 

ROS00650 

ROS00660 

R0SO0670 

RDS00680 

ROS00690 

ROS007CO 

R0S007I0 

R0S00720 

ROS00730 

R0S00740 

R0SO0750 

ROS00760 

RUS00770 

R0SO0780 

ROS00790 

R0S008O0 

ROSOOSiO 

R0S00820 

ROS00830 

ROS 00840 

ROS00850 

R0Sa0860 

BOS00870 

Rosooeeo 

R0S00890 

R0S00900 
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8 

n( I)=0(n+X(I) 



Rosoogio 


on TO 2 



ROS0092O 

9 

'^TRIAL=NTRIAL+1 



ROS00930 


PRUjT IO,F,{X{n ,1=1 ,N) 



ROS0O94O 

1.' 

FORMAK //1X,21H ♦^STARTING VALUeSX'*//7X 

,3HF *,F15.8//7X 

,3HX 

*, ROS0O95O 


lPF15.fi) 



ROS00960 


PRINT 11 



ROS0O97O 

1" 

FORMAT{//7X,3HV *) 



ROS0098D 


PRINT 12,((V(I,J),J=1,N),I = 1,N) 



ROS00990 

12 

FORMAT! 10X,4F15. 8/) 



ROS01000 

13 

FM=F 



ROSOIOJO 


IFCNSTAGE.LT.O) GO TO 19 



ROS01020 


DID =n(L)+Ya) 



ROS01030 


DU 14 I=1,N 



ROS01040 

14 

X( I)=0( 1) 



ROS 01050 


msl;xes=nsuxes+i 



ROS01060 


IF (NFFLUR.GT.O) GO TO 21 



ROSOIOTO 


E=ALPHA*E 



ROS 01 080 


GO TO 6 



ROS0109O 

15 

I F(LIMIT.eQ.O) GO TO 16 



ROSOllOD 


GO TO 17 



ROSOliiO' 

16 

IF(F.LE.FL) GO TO 17 



ROSOl 1 20 


GO TO 4 



RQSOilSO 

17 

IF(F.GE.FM) GO TO 13 



ROS 0.1 140 


GO TO 4 



ROS01150 

18 

IF(F.GE.FL) 60 TO 3 



ROSOildOi 


GO TO 4 



ROS01170 

19 

NSTAGe=NSTAGE+l 



Rosoiieo 


GO TO 6 



ROS01190; 

21 

YIL) =0.0 



Rosoizoo; 


L=L-i-l 



R0S012i0, 


IF(L.GT.N) GO TO 23 



R0S01220 

22 

F = S 



ROS0l23Ol 


NSuxes=o 



RUS01240 


NFELUR=0 



ROS 01 2 50 


GO TO 6 



R0S01260: 

23 

NSTAGE=NSTAGE+1 



ROSOl 270! 


CALL VECTOR! A, B,C,P,V,D,R,N) 



R0S012B0: 


TEMP1=R!1) 



ROS01290 ! 


TEMP2=0.0 



R0S013O0 


DO 24 1=1, N 



ROSOl 310 i 

24 

TEMP2=TEMP2+A!I ,2)?i'^2 



ROS01320 ; 


TFMP2=SQRT!TEMP2 ) /TEMPI 



R0S01330 : 


PRINT 25, NSTAGE, NTRIAL,F, TEMPI, TEMP2, (X!I), 1=1, N) 


R0S01340 ! 

25 

F0RMAT(//21H AT THE END OF STAGE ,I6,9H 

AND ,I6,9H 

TRIALS/ /7RDS01350 


1X,3HF *,F15.8//2X,8HMAG A1 # ,F15 .8 ,5X, 14HMAG A2/MAG Al*, 

F15. 

8//7X,R0S01360 ; 


23HX ’«',8F15.8) 



ROS01370 : 


PRINT 26 



ROS01380 ’ 

26 

FORMAT! //7X,3HV *) 



ROSOl 390 


PRINT 27,(!V(I,J),J=1,N),I=1,N) 



ROS01400 
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27 

FORMAlt 10X,4F15. 8/) 

ROS0141( 


IFINSTAGE.EQ.NS) GO TO 29 

ROSOL42( 


no 28 1=1 »N 

RUS0143£ 


n{ n =G.o 

R0S0I44C 


7{ n=n.o 

R0S0145{ 


Y{ I ) =0.0 

R0S0146( 


DO 2 8 J=1 ,N 

R0S0147C 


M I,J)=G.O 

R0S0148C 


3( I,J}=O.C 

R0S0149C 


C ( I, J)=0.C 

ROS0150<] 

28 

P{ I, JJ='J.C 

ROSOiSlO 


L = 1 

ROS01520 


GO TO 22 ^ 

R0S0i530 

29 

RC TURN 

R0SG1540 


FND 

ROS0155O 

C 


ROS01560 


SUBROUTINE VEC TOR { A ,B tC , P ,V , D, RtN) 

ROS0157O 


DIMENSION A(40,40) ,B(40,40) ,C(4U,40) ,P(40»4t)tV{40,40) 

ROS0158O 


DIMENSION D(40) ,R(4C) 

RaS01590 


DO 1 J=1,N 

ROS0160O 


DO 1 1=1 ,N 

Rosoi6:o 


DO 1 K=J,N 

ROS01620 

1 

A( I,J)=A{I,J) + D(K)*V(I,K) 

ROS0163O 


DO 2 1=1, N 

ROS01640 


B(I,1)=A(I,1) 

ROS0165O 

2 

R(1)=R(1)+B(I,1)**2 

R0S0166O 


R(1)=SQRT (RID) 

ROS0167O, 


DO 3 1=1, N 

ROS0168O 

3 

V(I,1)=B(I,1)/R(1) 

ROS0169O: 


DO 8 J=2,N 

ROS0170O: 


MJ=J-1 

R0S0171OI 


DO 4 K=1,MJ 

ROS0172O 


00 4 1=1, N 

R0S017 3Qi 

4 

CIJ ,K) =C( J,K) + A(I ,J)*va,K) 

ROS01740 


DO 6 L=1,MJ 

R0S01750: 


DO 6 1=1, N 

R0S01760! 


P(I ,J)=P(I,J)+C(J,L)*V(I,L) 

R0S01770; 

6 

B(I,J)=A(I,J)-P(I,J) 

R0S01780I 


DO 7 I=1,N 

R0S01790' 

7 

R{J)=RIJ1+B(I ,J)**2 

Rosoiaco; 


R{J)=SQRT (R{J)) 

Rosoiaioi 


DO 8 1=1, N 

ROSOl82Di 

8 

Vd ,J)=B(I ,J)/R{J) 

R0S0l830i 


RETURN 

ROS01840 1 


END 

ROS01850 j 

c 


ROS01860 ; 


FUNCTION FUNX0N{0,Z,NJ 

R0SO1870 i 

C 

VARIABLE XU) IS TEMPORARILY SUBSTITUTED BY VARIABLE Oil) 

R0S01880 ; 


DIMENSION 0(40) , Z(40) 

R0S01890 1 


FUNXON=100-* (0(2)-0(l)=t=0(l ) ) **2+ (l.-0( 1 ) ) **2 

ROS01900 



R0Sai91C 

R0SCi92C 

R0S0193C 

RUS0194fl 

ROS0195ti 
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APPENDIX D2 


GONEIDENQS INTERViiL CALODLi.TION FOR PiiRiiMETBR ESTIMi^TES 


The confidence intervals for the parameters were 
calculated from the ohservations as follows: 

1. The first order partial derivatives of the rate 
expression with respect to each of the parameters were obtained. 

2. The n-umcrical values of the derivatives were obtained 
by substituting the values of independent variables i.e. 
concentration and catalyst ratio and the estimated values of 
parameters. The values were written in the form of a nxm matrix. 
This nxm matrix is denoted by x where m being number of 
parameters and n being number of obseorvations. 

Let 

b^[RH]2 

® - bj+tjlM] 

be the rate expression. 

Differentiating with respect to b^, 

rsH i ^, 

- bj+b^TM] 


Similarly, the other two derivatives are 

1 2 


( /,R ) _ -h 

"'’2 (bj+bj [M]) 

i-.R [M] 

V 'T.-U / “ 




3 (b^+b^ [M]) 
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The matrix is of the form 


(-m.) 






X = 


^ dh^^n 


I dR N 
^ ^^2 n 


(^) 

^ 0b-,^n. 

j 


Where the subscript indicates the number of experimental 
observation. For example, in case of cumene a 14 x 3 matrix 
was obtained as 14 was the number of experimental observations. 

3, [x'^X] was obtained by inverting [X^X]. [X^X] ^ was 

a 3 X 3 matrix in all the cases. 

c. 


c = [X%]"^= 


11 


'21 


'31 


'12 


'22 


'32 


'13 


'23 


'33 


A typical [X^X]'"^ for the system cumene was, 


1 


18 

0.76112 

X 

10^5 

0.4426 X 10^^ 

,25 

0.46311 

X 

lo” 

- 0.26933 X 10^"^ 

,26 

-0.26933 

X 

10 

0.37569 X 10^^ 


' 0.12508 X 10' 

0.76112 X 10‘ 

0.44264 X 10' 

The parameter values for cumene obtained by Rosenbrock 
optimization procedure are as 


b 


2.2588 X 10' 
1.3744 X 10 

1,2048 X 10- 


10 
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The individual para,meter confidence intervals were 
ca.lculeted using the equations 


- h- 1 * Pi < h + h-f % 


’=2 - h - I \' f °22 £ p 2 < ^2 + h -| '''=22 


P’3 ^l- I ®yj_'''33 




2 
a 

'2 >^1 


where 


s- - estimated variance of y 


t a - value from 't* distribution corresponding to 
2 

to a significant level of a in a two sided 
test. 

s- was obtained from the exnression 


=^ 


RS3 

Y 


where RSS - residual sum of squares 
y - degrees of freedom 

= (number of obseivations - number of parameters). 



APPENDIX D3 
RESIDUAL ANALYSIS 

Table D,1; SrSTEM-CUMEUE 


Observed Rate 

7 

g mole/ sec x 10 

Predicted Rate 

7 

g mole/sec x 10' 

Hydrocarbon 
cone entration 
g mole/lit 

Residual 

X 10^ 

4.800 

5.030 

7.820 

-0.230 

4.130 

3,220 

6.260 

0. 910 

1.730 

1.810 

4.690 

-0.080 

1.140 

1.310 

3.990 

-0.170 

0.667 

0.805 

3.130 

-0.138 

4 . 380 

5.030 

7.820 

-0.650 

3.430 

3.220 

6.260 

0.210 

2.170 

1.810 

4. 690 

0.360 

1.220 

1.310 

3.990 

-0.090 

0.667 

0.805 

3.130 

-0.138 

3.370 

3.220 

6,260 

0.150 

1. 900 

1.810 

4. 690 

0.090 

1.500 

1.310 

3. 990 

0.190 

0.718 

0.805 

3.130 

-0.089 
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TABLE D-2’. SYSTEM-CYCLOHEIENE 


Observed Rate 

n 

g raole/sec x 10 

Predicted Rate, 

n 

g mole/sec xlO ' 

Hydrocarbon 

concentration 

g mole/lit 

Residual 

X 10^ 

4.90 

4.89 

7.88 

0.0119 

2.50 

2.76 

5.92 

-0.2590 

2.08 

1. 91 

4.93 

0.1670 

1.50 

1.23 

3.94 

0.2570 


TABLE D-3; SYSTEM-TETRALIE 


Observed' Rate,., 
g mole/ sec xlO 

Predicted Rate, 

7 

g mole/sec X 10 

Hydrocarbon 

concentration 

g mole/lit 

Residual 

X 10^ 

12.00 

12.00 

7.37 

0.0350 

7.84 

7.64 

5.89 

0.1960 

3.86 

4.30 

4.42 

-0.4440 

3.00 

2.99 

00 

• 

0.0081 




